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We show exactly with an SU (N) interacting model that even if the ambiguity associated with the 
placement of the chemical potential, [i, for a T = gapped system is removed by using the unique 
value /i(T — > 0), Luttinger's sum rule is violated. The failure stems from the non-existence of the 
Luttinger-Ward functional for a system in which the self-energy diverges. Since it is the existence 
of the Luttinger-Ward functional that is the basis for Luttinger's theorem which relates the charge 
density to sign changes of the single-particle Green function, no such theorem exists. Experimental 
data on the cuprates are presented which show a systematic deviation from the Luttinger count, 
implying a breakdown of the electron quasiparticle picture in strongly correlated electron matter. 



While the charge density remains fixed under 
renormalization from high (UV) to low (IR) energy, 
precisely what is carrying the charge can change 
drastically. For example, in QCD, free quarks at UV 
scales form bound states in the IR. The key signature 
of bound quark states is that the pole in the propa- 
gator is converted to a zero[l|, implying the fields in 
the UV-complete theory no longer propagate at low 
energy, hence, a breakdown of the elemental par- 
ticle picture. The conversion of poles to zeros of 
the single-particle Green function also obtains in su- 
perconductivity. In both QCD and superconductiv- 
ity, the new strongly coupled ground state in the IR 
lacks adiabatic continuity with the UV-state: free 
quarks (QCD) or free electrons (superconductivity). 
The question of how to compute the number of low- 
energy charged particle states is then problematic 
because what was a particle (pole) at high energy 
is no longer so at low energy. It is not surprising 
then that the only 'theorem', due to Luttinger Q, 
on the particle density in a fermionic system, makes 
no distinction between zeros and poles. The precise 
mathematical statement of Luttinger's theorem for 
spin-i fermions, 



= 2£)e(G(fc,w = 0)), 



(1) 



sums all momenta, k, where the Heaviside step func- 
tion, Q(G(k,uj = 0) is non-zero, with G(k,ui) the 
time-ordered single-particle Green function. The 
right-hand side of Eq. (Q} requires G(k,uj = 0) = 
or G(k,uj = 0) > 0. The latter obtains either from 
a pole or a zero of the single-particle Green func- 
tion. Hence, as far as the mathematics is concerned, 
poles and zeros of the Green function enter the par- 
ticle density on the same footing. We show that any 
statement of this kind in which zeros and poles are 
treated on equal footing is in error, hence the title 
of this paper. 



While poles of the single-particle Green function 
represent quasiparticles, zeros Q are quite different 
as they indicate the presence of a gap3- Equiva- 
lently, the self energy diverges, thereby represent- 
ing a breakdown of perturbation theory. As a 
result, purported non-perturbative proofs of Lut- 
tinger's theorem which assume gapless phases at 
the outset HOI 

bear no relevance to the validity of 
Eq. (JXJ) . Rather, such proofs are relevant only to 
the physical assertion that the volume of the Fermi 
surface is independent of the interactions — on some 
level a tautology, since all renormalizations from 
short-ranged repulsive interactions 043 are towards 
the Fermi surface. 

Naively, for gapped incompressible phases at T = 
0, the chemical potential can be placed anywhere 
in the gap with impunity. However, the placement 
of the chemical potential will change the energy at 
which G vanishes and hence affect the particle den- 
sity as pointed out first by RoschQ. Nonetheless, 
Farid[l(| has argued that the problem arising from 
this degree of freedom is entirely spurious because 
the chemical potential even at T = is unique, 
namely the limiting value of \i as T — > 0. For the 
case of the atomic limit of the SU(2) Hubbard model, 
this limiting procedure places the chemical poten- 
tial equally far from both bands, the particle-hole 
symmetric point, and Eq. (JTj) reduces exactly to 
n = 20(0) = 1, a result which holds beyond the 
atomic limit (Tlj . Farid's claim is interesting then 
because it would seem to establish a rigorous rela- 
tionship between a quantity which has no obvious 
physical import and a conserved one, the particle 



1 Trivial zeros that occur at the zone boundary as in the case 
of band insulators or symmetry broken states are irrelevant 
as they require two bands (see Appendix) whereas the zeros 
counted in Eq. JTJ are within a single band. 
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density. 

We show here exactly using an SU(N) general- 
ization of the atomic limit of the Hubbard model 
in which N flavors of "iso-spin" and n fermions re- 
side on each site that Farid's[l(J placement of the 
chemical potential does not salvage Eq. ([1]). The 
key result is quite simple. For this model, Eq. (fTJ) 
reduces to 



n = NQ(2n-N) 



(2) 



which is clearly not equal to the particle density un- 
less n e {0, N/2, N}. Namely, any partially filled 
band with TV odd leads to a violation of Eq. (fT]). 
That Eq. actually reduces to the correct result 
for the SU(2) case is entirely an accident because 
the O function only takes on values of 0, 1/2, or 1. 
The crux of the problem is that the Luttingcr-Ward 
(LW) functional strictly does not exist when zeros 
of the Green function are present. Since Eq. ([1]) 
relies explicitly on the construction of the LW func- 
tional and it does not exist for zeros of the Green 
function, Luttinger's theorem (Eq. (JTJ)) does not ex- 
ist. Our result applies to any system in which zeros 
occur, for example states of matter with a gap or 
a pseudogap, as in the cuprates. Experimental data 
on the cuprates is shown which support a systematic 
deviation from the Luttinger count. 

To illustrate the problem zeros pose for Eq. (JTJ), 
we consider for simplicity the SU(N) generalization, 



H = f (ni + n 2 H h n N f = f (hf 



(3) 



of the atomic limit of the Hubbard model. Since our 
key result hinges only on the existence of zeros, not 
on the details of a specific model, our conclusion is 
general. We have not included the site index here as 
it is superfluous to the many-body physics which is 
captured entirely by the N-flavors of "iso-spin" that 
five on a single site. Fig. (JTJ illustrates the level 
structure for N = 3. We define K = H(h) — fin and 
write the Green function, 



G a p(co) 



U 



2U 



2 



b ) + 



\ab 



(4) 
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FIG. 1. Schematic of the a) energy levels and b) spectral 
function for the Hamiltonian H — ^(h) 2 . For n = 2, 
e+ = H(n+1)-H(n) = 5*7/2, e_ = H(n)-H(n~l) = 
317/2, and as a consequence, /i = (e+ + £-)/2 = 2U. 



using the Kallen-Lehmann representation, where 
(a\c a \b) (b\c^\a) <a|c^|6> {b\c a \a} 



CO - Kh + K„ 



to - Kn + K h 



(5) 



in which the Green function is a sum of particle 
addition and removal parts. Here Z is the par- 
tition function and K a = (a\K\a). Since K a is 
completely determined by the occupancy number 
n a = J2 a ( a l c a c al a ) of the SU(/Y) orbital, we may 
write K a = K(n) for n a = n. Noting that there 
fN\ 

are I states of occupancy number n allows us to 



(6) 



simplify the Green function to 



N 



where 

Qa/3 = 



G a p(oo) = ^p(n)Ql 

n=0 



x +a p(n) 



x- a p(ri) 



lo - K(n + 1) + K{n) to - K(n) + K(n - 1) ' 



p(n) = i(^)e-^(»), ^ = Ep 
and the spectral weights are 



x +a p(n) = [ n ) E„ a =„ (a\c a cl\a) (7a) 



for particle addition and 



J (a\c^c a \a) (7b) 



for particle removal. 

To simplify the Green function, it suffices to calcu- 
late the spectral weights, x± a p. Note that in order 
for these matrix elements to be nonvanishing, one 
must have a = /3. In addition, the state a occur- 
ring in the summation expression for x+ a p (resp. 
X- a p) must be empty (full) at isospin a. There are 

^ n ^ lj^ sucn states, and so the final 

expressions for x± a p are 

x+ a p = S a p ( ' 1 ) / ( ] = 5 a p(l - f )(8a) 



and 



X-ap — fi<xf: 



N - 1 



8aP%. (8b) 



At T — 0, p(n) = 1 for some fixed n and p(n) = 
otherwise which effectively eliminates the sum over 
n in the Green function resulting in an expression of 
the form, G a p(oo) = 5 a[j Q n , where Q n (to) is Q2p(ui) 
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evaluated with Eq. ([5]). We now come to the issue of 
the chemical potential. According to Faridfioj. the 
leading /3 — > oo limit of the Green function is given 
by evaluating 

Q n (tu) + e-/ 3 ( Jf ("+ 1 )-*("))(Q»+ 1 (a;) - Q n (uj)) 

+ e^ K ^- K{n - 1 ^{Q n - 1 (io) - Q n (w)) 

+ o{e~ f){H{n+1) ~ 2H{n)+H(n ^ 1)) ). (9) 

The chemical potential is fixed by the relationship 



— = lim 

N P^oo 



Jr-i(-^ Giuj+[0) (1U) 



Combining the previous two expressions yields 

+ (1 _ 1L) e l3{K(n)-K{n-l)) 

+ 0(c~ l3{H{n+1) - 2H{n)+H(n - 1)) ), (11) 
which can be solved immediately to yield 
e + + e_ 



1 , N-n , _i, 

— log +o/3 1 

2(3 n 



(12) 



where e+(n) = H (n + 1) — H (n) and £-(n) = H(n) — 
H(n — 1). Consequently, the chemical potential is 
equidistant between the poles of Q n (ui). Equiva- 
lently, this choice for the chemical potential implies 
that K(n + 1) - K(n) > 0, K{n) - K(n - 1) < 0, 
&ndK(n+l)-K(n) = -{K{n)-K(n-l)). Hence 



G a p{u = 0) 



Sat 



K(n + 1) - K(n) 



2n-N 
A 



(13) 



Consequently, Luttinger's theorem for this system, 
if it is valid, is the statement that 



= ^9(G QQ (w = 0)) = A6(2n-A). (14) 



This expression clearly fails for any partial filling 
when A is odd. Also, by making the shift u — > 
uj + e(k)) in Eq. one sees that a band dis- 

persion arising from an infinitesimal hopping cannot 
change the argument of the 0-function in Eq. (fT4|) . 
Consider the atomic limit of the Hubbard model for 
A = 3 and two of the iso-spin levels with unit oc- 
cupancy, that is, n = 2 (see Fig. 1). This is the 
'half- filled' case. Eq. ([2]) clearly fails because &(x) 
can only take on values 1, or 1/2. Hence, no ex- 
pression of the form of Eq. ([2]) can ever yield the 
electron density when A is odd. At play here is the 
fact that particle-hole symmetry, which is present 
for A even, is strictly absent for A odd. 

Clearly if Eq. ^ fails, there must be an addi- 
tional term that contributes to the density. The ex- 
tra term is usually written as an integral involving 



the self energy. As (3 — > oo, G(oj) = G aa (uj) may be 
rewritten as 



G(w) = 



where 



uj + fi — e — E(w) 



E( w ) = "(^-") ( £ +- £ -) 2 



N 2 ui+^-eo 



and 



eo = ( f )e + +(l-f)e. 



(15) 

(16) 
(17) 

(18) 



The expression to be evaluated is 



h = -N lim - J2 GH^£HU iC (19) 



where the sum is over the fermionic Matsubara fre- 
quencies. The explicit calculation yields fioj 

h = ^ lim ((f)tanh/3(e+- M ) (20) 
+ (l - ^) tanh /3(e_ - fi) - tanh (3(e - fj,)) 
= 2 ( n s S n ( £ + - + - n) sgn(e_ - //) 
- Asgn(e -/•*))■ 

Because ll = e ~^ e+ and e_ < e+, we know that 
e_ < fx, e + > /z, and sgn(e — /i) = sgn(2n — A). 
Thus, 

J 2 = -(2n - A - Asgn(2n - A)) 



n - N<d(2n - A). 



(21) 



Combined with the previous result, I\ = A0(2n — 
A), we recover the full particle density, 



n = I\ 



(22) 



The failure of the LW identity, I2 — 0, rests en- 
tirely on the form of the self-energy in this problem, 
Eq. (fl7|) . Wc first note that E diverges at cj+li = e - 
Consequently, regardless of what is chosen for H , 
which in this case has been set to zero, E cannot be 
connected to any non-interacting problem as a re- 
sult of its divergence. Consider the LW functional, 
defined by 



SI[G] 
I[G = Go] 



duZSG 







(23a) 
(23b) 



which was used by LuttingerQ to show that the inte- 
grand of I2 is a total derivative. Because E diverges 
for some lo when G is the total Green function, it is 
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not possible to integrate the denning differential ex- 
pression in the neighborhood of the true Green func- 
tion, and therefore the LW functional does not ex- 
ist. Consequently, there is no Luttingcr theorem and 
Eq. ([I]) does not represent the density of a fermionic 
system because zeros of the Green function must be 
strictly excluded. While zeros are a robust math- 
ematical feature of a Green function, they do not 
represent a conserved quantity and do not have di- 
rect physical meaning. This can be seen from the 
fact that the zero line is sensitive to the placement 
of the chemical potential, and there is no observable 
consequence when the zero crosses /z. 

Even in the case where I2 = 0, and our system 
is gapped, the divergence of the self-energy is still 
present. Hence, no LW functional exists in this case 
as well. For N even, particle-hole symmetry is op- 
erative and it is this symmetry that results in a 
vanishing of I2 not any fact pertaining to the LW 
functional. For N odd, no such symmetry obtains. 
In gapless systems, Eq. (1) is still not generally 
valid. A less general formulation [TH, [l4[ which as- 
sumes the absence of zeros remains valid. That as- 
sumption, means that the interacting system must 
be perturbatively (adiabatically) connected to non- 
interacting fermions, which immediately rules out 
the Mott state which has a divergent self energy. 
In fact, the work by Hofava(l3j provides a promis- 
ing direction in which the robust features of a Fermi 
surface admit a K-theoretic formulation. 

The key implication the inapplicability of Eq. ([1]) 
portends for strongly interacting electron systems is 
that although the degrees of freedom that give rise 
to zeros undoubtedly contribute to the current, they 
have no bearing on the particle density. The particle 
density is determined by coherence (3E < e) while 
zeros arise from incoherence ( diverges [LL( sig- 
nifying that there is no particle to contribute to n). 
As a result, there exist charged degrees of freedom in 
strongly correlated electron matter which couple to 
the current but nonetheless cannot be given an inter- 
pretation in terms of elemental fields. Consequently, 
the particle density will be less than the total num- 
ber of degrees of freedom that couple to an exter- 
nal gauge field as demonstrated rece ntly[I3 when 
the upper Hubbard band is integrated out exactly. 
Note the breakdown of Eq. ([I]) has been demon- 
strated exactly in a model in which spin and charge 
are not separated, a purely atomic limit model in 
which there can be no difference between spin and 
charge velocities. 

Deviations from Eq. (Q} are expected then in 
experimental systems which are strongly correlated 
with either a hard gap or a pscudogap ( density 
of states vanishing at a single energy). Shown 
in Fig. ([2]) is a plot of the measured Fermi sur- 







+ 


s 




1 

+ 


+ / 

* 

/ 

/ 

/ 

/ 

/ 

/ 




X 


x=+V 






i 

s 

A 

t 








x Yang et al. (2011) 
+ He et al. (2011) 




1 1 



: I i L_ 

0.1 0.2 0.3 



x 

FIG. 2. Apparent doping zpg inferred from the Lut- 
tinger surface reconstruction as a function of the nominal 
doping x in LSCO and Bi-2212. 



face (area enclosed by the locus of k-points for 
which there is a maximum in the spectral func- 
tion) in Lai_ 2; Sr 2 ;Cu02[il| (+ plot symbol) and 
Bi 2 Sr 2 CaCu2 08+5 (x plotting symbol) as a func- 
tion of the nominal doping level in the pscudogap 
regime. A typical uncertainty in these experiments 
is ±0.02. Even when this uncertainty is considered, 
the deviation from the dashed line persists indicat- 
ing a breakdown of the Luttingcr count and hence 
a fundamental breakdown of the elemental particle 
picture in the cuprates. For the Hubbard model this 
systematic deviation has been seen previously pTj 
and stems from the fundamental fact that since the 
spectral weight is carried by two non-rigid bands, re- 
moving a single fc-state is not equivalent to removing 
a single electron. Another source of deviation from 
Eq. is the fact that as x nears optimal doping, 
the Fermi surface topology must change from scaling 
with x to 1 — x. Hence, there has to be a deviation 
xps = x - We hope this work serves to motivate a 
much more systematic study of deviations from Eq. 
(ID- 
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Appendix A 

The derivation of Luttingcr's theorem concerns 
only the eigenvalues of the propagator G on the 
Hilbert space of one-particle stateqj, not matrix ele- 
ments that happen to lie on the diagonal. The irrel- 
evance of the latter is underlined by the invariance 
of the entirety of the proof under the transforma- 
tion G i— > U GU^ , where U is an arbitrary unitary 
transformation on the Hilbert space. Consequently, 
methods of projection or integration of degrees of 
freedom that do not respect the eigenvalues of the 
propagator are incompatible with the familiar state- 
ment of the theorem, Eq. (1) in our manuscript. To 
illustrate, consider a model of two hybridized bands 
a, b. The Hamiltonian for this system is 

H = 22 (e ka a* k a k + e kb b* k b k + A k a* k b k + A* k b* k a k ) 

k 
k 

(Al) 

where 



tka Afe 

A fe e kb 



(A2) 



The propagator is then given by the 2x2 matrix 

1 

(A3) 



G(k,u) 



t*> - Uk 
1 



= u k 



oj-£k k 



where the diagonalization W. k = U k £ k Ul, U k lll = 1 
is given by 



U k 



Uk v k 



£k — 



£k+ 

£fe- 



(A4a) 



/ 2 i I A 1 2 - e ka + Cfcb 

e k ±= e k ± \/ e% + \A k \ , e k = , 

£fca — Cfcfi 



(k 



(A4b) 



Luttinger's theorem holds for this system: 

n = ^e(G ± (fc, W = 0)H^e(-e fc± ). (A5) 



fed 



fc± 



2 See, e.g., Dzyaloshinskii (2003). 

Note, Det[G] does not have any zeros; consequently 
any subsequent manipulations which result in zeros 
must be spurious. 

The aa component of G is 



G QQ (fc, oj) = 



\u k \ 



\Vk\ 



1 



u) - e k+ lo - £fc_ 



W — Cfca 



|A fc | 2 



(A6) 

and therefore the number density n a of electrons in 
band a is 



k 



(M 2 e(- £fc+ ) + |^|e(- efe _)). (A7) 



Suppose that the b k were integrated out. A naive 
substitution into the statement of Luttinger's theo- 
rem would yield 

n-i = ^e(G QQ (fc,w = 0)) = ^e(-£ fc+ £fc_/e fcfc ). 

k k 

(A8) 

The quantity n? is equal neither to n a nor to n. 
This failure stems from the fact that the Hamilto- 
nian is off-diagonal in the ab basis. In a general, 
non-diagonalizing basis, the propagator's matrix el- 
ements will be superpositions of poles and thus pos- 
sess zeros, none of which have anything to do with 
Luttinger's theorem. 



[1] G. 't Hooft, N. Phys. B 75, 461 (1974). 

[2] J. M. Luttinger, Phys. Rev. 119, 1153 (1960). 

[3] I. Dzyaloshinskii, Phys. Rev. B 68 85113 (2003); F. 

H. L. Essler and A. M. Tsvelik, Phys. Rev. B 71, 

195116 (2005). 
[4] M. Oshikawa, Phys. Rev. Lett. 84, 3370 (2000). 
[5] A. Praz, J. Feldman, H. Knorrer, and E. Trubowitz, 

Europhys. Lett. 72, 29 (2005). 



[6] J. Polchinski, arXiv:hep-th/9210046. 
[7] G. Benefatto and G. Gallavotti, J. Stat. Phys. 59, 
541 (1990). 

[8] R. Shankar, Rev. Mod. Phys. 66, 129 (1994). 
[9] A. Rosch, Eur. Phys. J. B 59, 495 (2007). 
[10] B. Farid, arXiv:0711.0952; arXiv:0711.3195; see also 

B. Farid and A. M. Tsvelik, arXiv:0909.2886. 
[11] T. D. Stanescu, P. W. Phillips and T. -P. Choy, 

Phys. Rev. B 75, 104503 (2007). 
[12] S. Hong and P. Phillips, arXiv:1110.0440. 



G 



[13] P. Hofava, Phys. Rev. Lett. 95, 016405 (2005). 
[14] G.E. Volovik, The Universe in a Helium Droplet 

(Clarendon Press, Oxford, 2003); M. A. Zubkov, 

arXiv: 1202.2524. 
[15] R.-H. He, et al, New. J. Phys. 13, 013031 (2011). 
[16] H. -B. Yang, J. D. Rameau, Z. -H. Pan, G. D. Gu, 

P. D. Johnson, H. Claus, D. G. Hinks, and T. E. 

Kidd, Phys. Rev. Lett. 107, 047003 (2011). 



[17] T. D. Stanescu and P. Phillips, Phys. Rev. B 69, 
245104/1-27 (2004); T. D. Stanescu and G. Kotliar, 
Phys. Rev. B 74, 125100 (2006); J. Kokalj and P. 
Prelovsek, Phys. Rev. B 75, 045111 (2007); B. S. 
Shastry, Phys. Rev. B 81, 045121 (2010); for the pe- 
riodic Anderson model, see Z. Gulacsi and D. Voll- 
hardt, Phys. Rev. B 72, 075130 (2005). 



